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Abstract 

The semiclassical quantization conditions for all partial waves are derived for bound states of 
two interacting anyons in the presence of a uniform background magnetic field. Singular Aharonov- 
Bohm-type interactions between the anyons are dealt with by the modified WKB method of 
Friedrich and Trost. For s-wave bound state problems in which the choice of the boundary condi- 
tion at short distance gives rise to an additional ambiguity, a suitable generalization of the latter 
method is required to develop a consistent WKB approach. We here show how the related self- 
adjoint extension parameter affects the semiclassical quantization condition for energy levels. For 
some simple cases admitting exact answers, we verify that our semiclassical formulas in fact provide 
highly accurate results over a broad quantum number range. 



* Electronic address: hurjin2@snu.ac.kr 



^Electronic address: cklec@phya.snu.ac.ki 



I. INTRODUCTION 



In two-dimensional space where the rotation group (SO (2)) is abelian, particles obeying 
fractional statistics, i.e., anyons, can exist [|l|, |^, |[]. They can be given a concrete mathe- 
matical representation in flux-charge composites, since 'charged' bosonic (or fermionic) parti- 
cles with non-zero magnetic flux attached behave like anyons thanks to the Aharonov-Bohm 
type interference ||. Among the physical realizations of anyonic objects, the quasiparticle 
or quasihole in the Laughlin state of the fractional quantum Hall system corresponds to a 
particularly interesting case M. In this regard, a challenging problem is to settle from the 
first principles what sort of quantum states would be allowed for anyon matter with some 
realistic interactions (and possibly a uniform background magnetic field). 

Choosing bosonic description, quantum dynamics of a system of anyons will be governed 
by the Hamiltonian of the form || 
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H = J2 — p n -aa n -eA e *(f n ) + - £ V(\r n -f m \), (1.1 
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where r n = (x^,x^), v4 ex (r) = — (— x 2 ,x l ) with B > is the vector potential appropriate 
to a uniform background magnetic field, and 
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The parameter a characterizes the type of anyons and without loss of generality a may be 
restricted to the interval (-1, 1]. For a = (a = 1), we have physical bosons (fermions). This 
defines a Galilean-invariant system if the background magnetic field, included for generality, 
is turned off (i.e., for 5 = 0). In the presence of a generic two-body interaction potential 
V(\f n — f m \) together with the Aharonov-Bohm interactions represented by the point vortex 
potentials ( |1.2|) , the above Hamiltonian leads to a nontrivial energy eigenvalue equation even 
for the two-body case — the starting point of a systematic n-body study. [See Refs. [|3], [|, 
0, §, ^ where some related studies are made, usually without a two-body potential term]. 
For the two-anyon system the time-independent Schrodinger equation, after separating the 
center-of-mass dynamics (just the problem of a particle with mass M = 2m and charge 2e 
moving in the background magnetic field), reduces to the following equation which has only 
the relative position r = r 2 — r*i = (x, y) as independent variables: 
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+ V{r)\i}{r) = Ei<(r), (1.3) 
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where fi = —m and a l (r) = ^ — ■ [We set % = 1 in this paper]. The relative dynamics 
is effectively that of a single particle moving in the presence of a uniform magnetic field, a 
point vortex at the origin, and a certain radial potential V(r). 

With a nontrivial potential V(r) an exact analysis of (|1.3|) is usually not possible and 
hence suitable approximation methods may be sought. In this paper, we shall study quan- 
tum bound states of two anyons through the semiclassical or WKB analysis of the radial 
Schrodinger equations for partial wave amplitudes ipi(r) (I = 0, ±2, ±4, • • •) which are de- 
rived from ( |1.3|) . To account for the effects due to the singular point vortex potential (as 
well as the centrifugal potential term) within the WKB method, the conventional Langer 



modification [10| of the potential is not adequate; but, the modified method of Friedrich and 



Trost [ I I ] can be used in an effective way. Based on the latter method, we have obtained 



the semiclassical quantization condition (containing a nonintegral Maslov index fT2|j ) 

jH drp^E-V eff {r)) = (n + ^ + ~|J - a\ - ^) ^ ( n = 0, 1, 2, • • •) (1.4) 

where r\ = (I — a) 2 — - (assumed to be positive here), V e g{r) is the effective one- dimensional 
potential 

Mr)= V (r) _^ (I - a) + |. B V + ^ 5 , (1 .5) 

and ri and r2 refer to the related classical turning points. [The condition ( |1.4|) may be used 
when the WKB wave function in the region r > r2 contains a decreasing exponential only]. 
Efficacy of this formula becomes evident once one sees how its predictions compare against 
the exact values available for some simple cases. 

Actually, for the case of the s-wave amplitude ipi=o(r), there exists an additional compli- 
cation involving the choice of boundary condition at r = (i.e., at the point of two-particle 
coincidence). We know from the theory of self-adjoint extension that there exist a one- 
parameter family of acceptable boundary conditions fl3| , including the so-called hard-core 
boundary condition || as a special case. [Note that there is no a priori reason to choose 
specifically the hard-core boundary condition — a real system under study should deter- 
mine the relevant boundary condition eventually]. In accordance with this theory, one can 
represent the s-wave amplitude ipi=o(r) for small r by the form 

V>i=o(r) oc I J\ a \(kr)+ tan0 (^\ J^ a \(kr)\ , (1.6) 



where k = J2p{E — V(0) — 9 (the self-adjoint extension parameter) is a real dimen- 

sionless number, and p a reference scale introduced for convenience. The hard-core boundary 
condition is associated with a special choice, 8 = 0. Needless to say, this boundary-condition- 
dependent effect should be taken into account in the WKB analysis of the s-wave amplitude. 
The Langer modification method is simply not applicable here; but, the method of Friedrich 
and Trost has a natural generalization for this problem. The resulting s-wave semiclassical 
quantization condition we have obtained for — < \a\ < 1 has the form 

jT dr^2p(E-V eff (r)) = (n + ^ + ^\a\-y^-^j7t + e(E), (n = 0, 1, 2, • ■ •) 

(1.7) 

where 0(E), a function of energy E (through k = J2p~^E — V(0) — j^ot) ), is related to the 

above self-adjoint extension parameter by 

^ sin7r|a| tan# (k/p) 2 ^ 

1 + cos7r|a| tan# (k/p) 2 \ a \ v ' 

For I a I < -, we need some extra consideration and this case is covered by our another 
2 



formula, given in ( |3.12|) . 



It might be of some interest to look at our problem also from the viewpoint of field 
theory. In the field-theoretic description, one can describe (nonrelativistic) anyons by using 
bosonic Schrodinger fields i/j(r,t), ip^{f,t) coupled to an abelian Chern-Simons gauge field 
a M (r, t) (/i = 0, 1, 2) |14|| . For the above system specifically, one may consider the Lagrangian 
density 

--J d 2 f'^\f,t)^\f',t)V{\r-r'\)-ilj{f',t)^{f,t), (1.9) 

where D t = d t + iqao, D = V — iqa — ieA cx , and the anyon parameter a should be identified 

q 2 1 

with . Here the contact interaction term, — X^ib ip ipip, which becomes necessary to 

2iik 2 

ensure the renormalizability of the theory, is responsible for the s-wave boundary condition 
ambiguity mentioned above |15|]; for the precise correspondence between the method of self- 



adjoint extension in quantum mechanics and the regularization/renormailization procedure 



in field theory, see Refs. JT^, 0, PH . Thus, with a suitable transcription made from the 
self-adjoint extension parameter to the renormalized contact coupling \r (as considered, for 
instance, in Ref. fl?l), our WKB analysis provides appropriate results for the two-particle 



bound states of this nonrelativistic field theory system also. But, in this work, such field 
theoretic language will not be used. 

This paper is organized as follows. In Sec. 2 we use the modified WKB method to find 
approximate non-s-wave bound states of the Schrodinger equation ( |1.3j ) and especially derive 
the semiclassical energy level formula ( |1.4j ). Then, for some simple cases (e.g., for a circular 
billiard), we check our WKB-based predictions against the exact results. In Sec. 3 the case 
of s-wave bound states is studied within the framework of the modified WKB method, with 
special attention given to the dependence on the self-adjoint extension parameter (or contact 
coupling). Section 4 contains a summary and discussions of our work. 



II. NON-S-WAVE SEMICLASSICAL BOUND STATES 



The relative dynamics of the two anyon system is governed by the Schrodinger equa- 
tion ( |1.3| ). This equation is analyzed most conveniently in polar coordinates, as the vector 
potentials in our problem are equal to 



a(r) = -9 
r 



A cx (r) = -r§. 
v ; 2 



(2.1) 



Now, writing ip(f) = il)i(r)e lW (with I restricted to even integer values in our bosonic de- 
scription of anyons) in ( |1.3|) , one obtains the following radial equation for the partial wave 
amplitude i^i{r): 
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where V(r) = V(r) (/ — a) H r 2 . Introducing the function Ri(r) by 

4:[I 32/1 

if>i(r) = —=Ri(r), 



( f2.2| ) can be further recast into the form 

f 1 d 2 
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+ V eff (r)\R l (r) = ER l (r) 



with the effective one-dimensional potential (given already in (|1.5Q ) 



(2.2) 



(2.3) 



(2.4) 



Ve ff (r) = V(r) + 
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{r, = {l-a) 2 --). 



(2.5) 



The eigenvalue problem (|2.4|) for I = ±2, ±4, • ■ • (and hence 77 > 0, assuming \a\ < 1) will be 
studied by the semiclassical method in this section. The I = case is considered separately 
in the next section. 

A naive application of the WKB method with the differential equation ( |2.4| ) does not 
yield very satisfactory results (especially for relatively small I), the reason being that (i) the 



radial coordinate r runs from to oo (instead of — oo to oo) and (ii) the potential contains a 

2 fir 2 

the correct small-r behavior. This defect gets significantly reduced if one introduces, within 



singular term at the origin, - — -. In fact, the naive WKB wave function does not even have 



the usual WKB approach, the so-called Langer modification of the potential [[[0 , effected 
through replacing the potential (|2.5| ) by 

V + z 

"« l, ' , = l in+ 2^" 



Vt ff {r) = V{r) + (2.6) 



But, to go beyond this simple Langer modification scheme, it becomes necessary to incor- 
porate the correct phase loss due to reflection at a classical turning point in the WKB wave 
function corresponding to the classically allowed region. This leads to the modified WKB 



method, as described recently by Friedrich and Trost JTT|. In this paper, we shall use the 
latter method for our semiclassical discussion and see how the resulting predictions compare 
against exact results (and also those obtained with the Langer modification) for some special 
cases. 

We may here suppose that the potential V e ff(r), as given by ( |2.5| ) (with r] > 0), has a 

typical shape shown in Fig. |I[ Also it may be assumed that the potential V e ff(r) near the 

~ ?7 ~ eB 
origin can be approximated by V(0) H -, with V(0) = V(0) (/ — a). Then the exact 

solution of (|2.4|), which is regular at the origin, should have the small-r behavior as given 



by 

R^(r) oc VkrJ v (kr), (2.7) 



where k = \j2fx(E — V"(0)), and J v is the Bessel function of order u = Jr] + j — \l — a\. We 
are here interested in finding the WKB bound states in the situation where there are two 
classical turning points (see Fig. [I]), ri and r 2 . The WKB wave function in the classically 
allowed region r\ < r < r 2 may be written as 

RT KB (r) oc -jL= cos { f dr'p(r') + <j)) , (2.8) 
./p(r) ' 
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FIG. 1: The shape of our potential V e g(r) with r\ > 0. 



where p(r) = ^2p,(E — V e g(r)) = ^jk 2 — \ ~ 2£t(^( r ) ~ ^(0)), and the phase is yet un- 
determined. Note that 20 can be identified as the phase loss of waves due to reflection by 
the barrier at r < r\. In the conventional WKB method, one then uses the famous connec- 
tion formula to combine the oscillating wave function in Q2.8Q with suitable monotonically 
decreasing wave functions in the classically forbidden regions (i.e., < r < r\ and r > r 2 ). 
But, by the reasons mentioned above, a different strategy must be adopted in our case in 



dealing with the WKB wave function for small values of r. Here, following Ref. ||11|| , we will 
simply fix the phase <fi in ( |2.8| ) in such a way that the asymptotic expansion of the expression 

(i.e., the correct expression 



Q2.8|) for relatively large kr may match that following from 

when r is not too large), say, up to terms of order — — — . On the other hand, in the vicinity 

(kr) 2 

of the second turning point r = r 2 , the normal WKB method with the standard connection 
formula may well be used. 

Using ( |2.7|) , the exact solution of the Schrodinger equation near the origin should have 
the asymptotic expansion 
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(kr) 



(2.9) 



For the asymptotic expansion of the WKB wave function (|2.8|) , we may use the small- 

r approximation of V e ff(r), i.e., V e ff(r) ~ V(0) H in p(r) (and correspondingly the 

2fir 2 
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turning point value r\ — -^—), and then the integral / dr'pir') can be performed explicitly. 

k Jn 

The result is the following asymptotic behavior: 



R7 KB (r) ~ I 1 



Tj{fj 



\(kr) 



cos [kr — c + 



2 At 



Tj ( 1 

sin [kr — c + 0] + O 



(kr) 



(2.10) 
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where c = ■ Comparing ( |2.10| ) with (j2.9|), we see that the two expressions agree up to 



terms of order — — 777 only if we choose the phase <j) as 



(kr)' 



c- [u + 



1\ 7T 



2/ 2 



1\ 7T 

27 2' 



(2.11) 



[In contrast, we note that the standard WKB method with the Langer modification pro- 



duces the right argument of sine and cosine in (2TT0), but wrong coefficients for the terms 

1 1 
proportional to — and — — — in (|2.10|)]. 



kr (kr) 2 

The phase <fi in our WKB wave function (|2.8| ) has been determined now. On the other 
hand, the WKB wave function obtained by applying the standard connection formula at the 
second turning point r = r2 has the form 



Rr KB (r) oc 



p(r) 



cos 



ra 7T 

drp(r) + — 



(2.12) 



The two functions ( |2.8| ) and ( |2.12| ) must of course be the same. From this follows the 
quantization condition for energy levels 



ra 7r 
dr'p(r') = nn — (f) -\ — 
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(n = 0,l,2,---) 



(2.13) 



This is our formula ( |1.4|) . For the sake of comparison, we note that the WKB quantization 
condition obtained with the help of the Langer modification reads 



[•T2 




/ dr. 




>ri \ 
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it. 



(n = 0,1,2, 



(2-14) 



We also remark that, for an effective potential having a qualitatively different form from 
that shown in Fig. [I], the WKB wave function beyond the point r = r2 might not be given by 



7T 



a decreasing exponential only. In such case the phase — in the argument of cosine in (|2.12|) 



should be replaced by an appropriate different value, and our formula for the semi classical 
energy levels need to be adjusted accordingly. 

We will now consider some special cases to see how well the predictions based on our 
semiclassical formula ( |2.13| ) (and fl2.14| ) also for comparison) fare against the exact results. 
Suppose we have V(r) = 0, but B ^ 0, i.e., the case of noninteracting anyons in a uniform 
magnetic field 0, 0]. In this case the solution to the radial equation (|2"^|), which is regular 
at the origin, can be expressed in terms of the confluent hypergeometric function fT8H , 



b 1 



Ri(r) = (^r 2 ) " e-^ 2 / 8 F(a, b, ^r 2 ), (2.15) 
where a = -[1 + \l — a\ + (Z — a)} and b = 1 + \l — a\. The large-x asymptotics of 

~" r(^) 

the confluent hypergeometric function F(a,b,x) is well-known: F(a,b,x) ~ e x x a ~ b if 

T(a) 

a/0,-1, —2, ■ ■ •, while F(a, 6, x) for a = 0, —1, —2, • • ■ reduces to a polynomial (called the 
associate Laguerre polynomial). Hence, to obtain a normalizable solution from ( |2.15| ), the 
energy E must be restricted to the values 

(n + i + - a\ - - a)) u e , (n = 0, 1, 2, • • •) (2.16) 
eB 

where lu c = — is the classical cyclotron frequency. These are exact energy levels, exhibiting 
very different /-dependences for / > and for Z < || On the other hand, within our 

semiclassical approach, the energy levels follow immediately once we evaluate the integral 

eB e 2 B 2 

in the left hand side of Q2.13| ). With V(r) = (Z — a) H r 2 , the given integral is 

Afx 32/i 

readily evaluated to yield the condition of the form 



E 



= n + l + K^V + l-Jv)}*. (n = 0,l,2,...) (2.17) 

If we solve this equation for E, the result is precisely fl2.16|) ; in this case, our semiclassical 
formula for energy levels yields the exact results. Incidentally, for this special case, we note 
that the exact results are reproduced even if we use the WKB formula with the Langer 
modification, i.e., (|2.14|) . 

As a particularly simple example with nonzero interaction, let us now consider a circular 
billiard potential (for general B > 0), i.e., 

' , r < R 

V(r) = { (2.18) 
oo , r > R. 



Veff{r) 



E = E r 




E < E r 



FIG. 2: The shape of V e jj{r) for a circular billiard of radius R. 



[This type of potential was previously considered in connection with the virial coefficient 
calculation of the anyon gas PI|]. The corresponding exact solution to Q2.4|) is then the 
expression ( |2.15| ), subject to the boundary condition R[(r = R) = 0, that is, 

F (a,b,z= ^-R 2 ) = 0, (2.19) 

when a and b here represent the same quantities as defined above. Given information on 
the roots of (|2~T9|) , one can determine the exact energy eigenvalues. On the other hand, for 
the semiclassical energy levels, we cannot use the formula (|2.13|) blindly — for the effective 
potential as shown in Fig. 0, the expression ( 2.12 ) is not appropriate. That is, instead of 
( [2.12j ), we must use the form 



R, 



WKB ( 



r) oc 



p(r) 



cos 



r-2 



dr'p(r') + 5 



[2.20) 



7T 



with the phase 5 ( ^ — ) calculated to suit the present case. This in turn has the consequence 



that we must add [5 — —J to the right hand side of the quantization condition ( p,13|) (and 

of ( |2.14| ) also). Here we have two different situations depending on E > E CT or E < E cr , 

eB e 2 B 2 ri 

where E CT = (/ — a) H R 2 H -. For E > E CT , the second turning point is at 

4/x 32/i 2jj,R 2 

7T 

r 2 = R (i.e., at the infinite wall) and, obviously, 5 = — is the correct choice. For E < E cr , 
however, we have r 2 < R (see Fig. 0) and may consider another WKB wave function in the 
classically forbidden region r 2 < r < R, 



R7 KB (r) oc 



Ae Jr 2 



SI dr'K{r') + B J r r2 dr>K(r>) ^ ^ _ / 2 ^V eff (r) - E) (2.21) 



10 



TABLE I: Energy levels (E n ) for a circular billiard of radius R with nonzero B. The parameters 

1 1 

were chosen as I = 2, a = — , fj, = — , R = 1, eB = 2. 



n 


Langer 


Ours 


Exact 





18.4719 


18.7389 


18.7843 


1 


57.9961 


58.2515 


58.2664 


2 


117.2249 


117.4775 


117.4850 


3 


196.1862 


196.4378 


196.4425 


4 


294.8844 


295.1355 


295.1385 


5 


413.3210 


413.5717 


413.5745 


6 


551.4962 


551.7468 


551.7484 


7 


709.4104 


709.6608 


709.6621 



with the coefficients A and B chosen such that i?^ VKB (r = R) — may hold. Then the 
phase 5 in ( [2.20|) should be determined by using the (standard) WKB connection formula 

7T 1/2 f R dr'ntr') \ 

with (|2T2Tl) ; this gives the value 5 = - + cot -1 I 2e Jr 2 k ; 1 . 

In the presence of the billiard potential fl2.18|) , we have given above all the elements 



from which the exact and semi classical energy levels can be found, numerically at least. By 
studying them, we have observed that the semiclassical predictions, especially the ones based 
on our method (compared to that based on the Langer modification of the potential), give 
generically highly accurate results for a broad range of quantum numbers. This should be 
evident from looking at Table |, where the exact levels and the semiclassical predictions - 
both ours and that based on the Langer modification — are compared for some representative 
choice of parameters (for which we have r2 = R for all n — 0, 1, 2, • ■ •). [With R = 5 (but 
the same values for the other parameters) we have r 2 < R for small n; still, we have found 
the excellent overall agreement between the semiclassical and exact energy levels]. 

III. S-WAVE SEMICLASSICAL BOUND STATES 

In this section we specialize to the s-wave bound states of two anyons. The relative 
wave function then depends on the coordinate r only, i.e., ip(r) = —=R (r), and from the 



11 



Schrodinger equation the function Ro(r) should satisfy the differential equation (see ( |2.4j )) 

f i£ a 2 - i ~ 1 

+ ^-~± + V(r)\R (r) = ER (r), (3.1) 



[ 2/i dr 2 2/ir 2 } 
eB e 2 B 2 

where V(r) = V(r) H a H r 2 (the / = form of the same function introduced in 

4/x 32/i 

the previous section). For any potential V(r) which is regular at r = 0, one may infer the 

small-r behavior of the solution Ro(r) by studying the solution of the simpler second-order 

differential equation, obtained from ( |3.1|) with the replacement of the potential V(r) by 

eB 

the constant V(0) = V(0) + - — a. In fact, it was an analogous reasoning that picked the 
form Q2.7D as a unique allowed choice for 1^0. For s-wave states, on the other hand, the 



form ( |2.7| ), i.e., Rq x (t) oc \fkrJ\ a \(kr) does not represent the most general small-r behavior 
allowed for the exact solution of fl3.1|) ; here, one finds a normalizable wave function even if 
Ro(r) contains a piece corresponding to another independent solutioin of the just-mentioned 
second-order differential equation, that is, the form y/krN\ a \(kr) (N\ a \ is the Neumann func- 
tion). Therefore, for the small-r behavior of the exact solution of (|3.1| ), one may well consider 
the general form 

Rl x (r) oc Vk^ {cosG J\ a \(kr) + sinQ N H (kr)} , (k = ^2fi(E - V(Q)j\ (3.2) 

with an arbitrary r-independent angle G. 

In quantum mechanics with a singular potential, it is the theory of self-adjoint extension 
of the related Hamiltonian that determines what boundary conditions may be allowed. [A 
particularly instructive example here is that of the two-dimensional 5-function potential 



problem, considered in Ref. [19]]. From this theory follows the one-parameter family of 
boundary conditions for our s-wave function ipi =0 (r) = —=R (r) [|T3||, 



lim {rW^(r) - tanO g^'"' ffi±H)_^. [rW ^ (r)] j = , (3.3) 

where 9, a dimensionless real number, is the self-adjoint extension parameter, and p just a 
(conveniently introduced) reference scale. It is this general boundary condition that leads 
to the small-r behavior shown in ( |1.6| ). One can now translate the expression ( |1.6| ) into 
the form ( |3.2| ) with the help of the identity J^\ a \(kr) = J\ a \(kr) cos7r|a| — N\ a \(kr) sin7r|o;|. 
Between the constant G in (|3.2|) and the self-adjoint extension parameter 9, we then find 
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the connection given in (|1.8|) ; this also tells us that depends on k (and hence on energy 
E). 

What would be the correct way of implementing the specific small-r behavior shown in 
( p.2|) with semiclassical bound states? Clearly, applying the standard WKB method with 
the already-mentioned Langer modification trick would be unjustified here. But we can still 
make use of the Friedrich-Trost approach — use the WKB wave function in the classically 
allowed region, i.e., 

1 



#™(r) oc cos \ I dr'p{r') + <f>\ , (3-4) 

lp(r) 



p(r) = JME - V eff (r)) ; V eff (r) = V(r) + 



2\xr 2 J 

with the phase <fi chosen such that this WKB wave function may lead to an asymptotic form 
(for relatively large kr) in agreement with that following now from the supposedly correct 
expression in ( p.2|) . The asymptotic expansion of (|3.2|) is easily found: 

' sin !kr-(\a\ + l)^--e\+o( —!— J , (3.5) 



2kr I V 2 J 2 J \(kr 

where r] = a 2 — —. As regards the WKB wave function (|3.4j) , a new problem arises for 

i] < (i.e., if V e ff(r) contains an attractive singular term at the origin) and it thus becomes 

1 . 1 

necessary to discuss the case with |a| < - separately from that with \a\ > -. See below on 

this. 

If | a | is lager than - (i.e., for 77 > 0), the present s-wave WKB wave function will have 



qualitatively the same structure as the non-s-wave WKB function considered in (|2.8|) , and 
hence its asymptotic behavior by the form ( |2. 10|) with the substitution rj = a 2 — — now. 



4 



Then, as we demand that it match the asymptotic behavior in ( |3.5|) , the phase in the 



WKB wave function ( |3.4| ) will be fixed to have the value 



f f 2 - H |a| + i)i - e - ( |a|a 5^ ,3 '" 1 



[Note that, with B set to zero, this is nothing but the formula ( |2.11D . Also we have invoked 



the continuity of our formula to include the case \a\ = — here]. If |a| < - (i.e., for 77 < 0), on 
the other hand, we find the behavior V e g(r) — > —00 as r — > and hence no classical turning 
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Veffir) 




FIG. 3: The shape of V eff {r) with 77 < 0. 
point near the origin. (See Fig. |3|). Thus, r\ in the expression ( |3.4| ) is ambiguous. For this 

rr 

case, one might be inclined to choose n = 0; but, then, the integral / dr'pOr') becomes 

Jo 

ill-defined. This difficulty is resolved if one chooses the point r\ to be not the origin but a 



certain (small) positive number [20]. The asymptotic behavior of the WKB wave function 
(p.4|) can then be described by the form (|2.10|) again only if the constant c there is now 
taken to be 



\J\f]\ + (kri) 2 — J]rj\\n 



\i\ , \/\<n\ 1, 2 1 



f + W^ + "^ • {^^a'- l <0). (3.7) 



Comparing such asymptotic expressioin with the behavior ( |3.5| ), one immediately sees that 
the phase in this case should be chosen 

«l + i)i-e, (H<i). (3.8) 

The WKB wave function ( |3.4|) for \a\ < —, with this phase choice will have (in the region 
beyond the point ri) only a very mild dependence on the cutoff r±. Actually, we can do 
away with introducing this ad hoc cutoff by changing our WKB wave function slightly. Let 
V e ff(r) denote the small-r approximation of V e ff(r), i.e., 

V eff (r) = V(0) - (3.9) 

Then we would not sacrifice none of the desired properties of the WKB wave function by 
changing the argument of cosine in (|3.4|) to the expression 



f Q dr> { ^(E-V eff (r')) - p^E-V eff {r'))^ + jT dr' ^(E - V eff (r')) + 
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dr' J2h(E - V eff (r')) - ME - V eff (r')) + VM + (At) 



— J \t)\ In 




Oil 



71 

2) 2 



-e, 



(3.10) 



has been used. With |a| < -, this manifestly cutoff independent WKB wave 



where 

function shall be used in our discussions below. 

In the presence of the second classical turning point r = r 2 , the WKB wave function 
found above can be used to obtain the s-wave semiclassical energy levels. This part of 
argument is entirely parallel to that in the previous section, and hence we may state the 
result only. Assuming that the potential is smooth near the second turning point and 
the WKB wave function in the region r > r 2 contains a decreasing exponential only, the 
quantization condition for energy levels reads, for \a\ > - (i.e., r] > 0) 



2 drp^E-V eB {r)) 

and, for |a| < - (i.e., — - < rj < 0) 
fj dr [sjME^V^) - 



1 1 

n H h - 

.22 



vr + e, 



(3.11) 



— v I *7 I m 



1 + 




;jff (r))J + y/\rj\ + (kr 2 y 



n 



\a\ ix 



e. 



(3.12) 



kr 2 I V 2 2' 

Here it should be remarked that G itself is a function of E, being related to the energy 
variable (and the self-adjoint extension parameter 9) by ( |1.8| ), It is through the presence of 
this function G in ( |3.11| ) and ( |3.12 ) that the ^-dependence enters the semiclassical energy 
levels. 

To check the accuracy of the above s-wave quantization formulas, we may here consider 

1 

the case V(r) = 0, B ^ (and with \a\ > — for simplicity). Then the exact solution to 
J]), which satisfies the 6*-dependent boundary condition (|1.6|) at the origin, is 



Ro(r) 



Z 2 4 e 2F(a,b,z) 



2\a\ 



+tan "U U 



/eB\ H r(l + \a\) 



i _ eB 2 

where z = — — r , 
4 



r(i-H) 

2fiE 



2 (1 + |a| + a) - eB 



z~2e-2F(l + a-b,2-b,z), (3.13) 
and b = 1 + lotl. For this function to be 



normalizable, we must further demand from the asymptotic consideration that 



r(a) 



2\ 2|Q| (eB\ H r(l + |a|) r(2 - b) 



tan 9 



\ 4 J r(l- |a|)r(l + o-6) 



0. 



(3.14) 
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TABLE II: Energy levels (E n ) for the s-wave of the eigenstates of V = 0, but B / 0. The 

2 1 7T 

parameters were chosen as a = — , u = — , eB = 2, 9 = — . 

3' P 2' ' 4 



n 


Ours 


Exact 


-1 


1.6667 


1.6477 





3.3468 


3.3440 


1 


5.2337 


5.2329 


2 


7.1796 


7.1793 


3 


9.1480 


9.1478 


4 


11.1270 


11.1269 


5 


13.1120 


13.1120 


6 


15.1007 


15.1007 



Solving this equation for E will give exact energy levels. We may then compare the semi- 
classical energy levels obtained with the help of ( |3.11 ) against these exact values. [With 
9^0, the semiclassical formula does not produce exact energy levels]. In Table [HJ, such 
comparison is made for some representative choice of parameters. We find the high accuracy 
of our semiclassical predictions, even for small n, very impressive. 



IV. SUMMARY AND DISCUSSIONS 



In this work we have developed the semiclassical theory of two-anyon bound states, with 
careful consideration given on the treatment of singular Aharonov-Bohm-type interactions 
between anyons. The modified WKB method of Friedrich and Trost has proved to be 
particularly effective for this problem. We have also clarified the role of the self-adjoint 
extension parameter for s-wave bound states within this semiclassical approach. For some 
simple cases, we have been able to confirm that our semiclassical formulas provide highly 
accurate energy levels over a broad quantum number range. We expect this to be the case 
with more general classes of interaction potentials. It should also be possible to extend 
this semiclassical theory to the case of two-anyon scattering states (in the absence of a 
background magnetic field). 
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Results of the present semiclassical theory may be applied to study certain features of 
the anyon gas. In the high-temperature low-density limit, for instance, one usually resorts 
to the virial expansion to study various thermodynamic properties. In the case of the 'free' 
anyon gas, the second virial coefficient has been calculated in Refs. [F2T|, |2"2], B3|; especially, in 



Ref. J2"3"|, the effect of the self-adjoint extension parameter on the virial coefficient has been 



considered also. Based on the semiclassical understanding of two-anyon bound/scattering 
states, one may extend this discussion to the case of the anyon gas with some nontrivial 
2-body potential, and see what new features the presence of such 2-body interaction can 
give rise to. Some of these issues are under investigation. 

Also note that more general kinds of anyons, other than the ones we discussed here, are 
possible. Anyons obeying so-called matrix (or mutual) statistics [f24], [25[ may prove to be 
relevant in the discussion of multi- layered quantum Hall effect, and, theoretically, particles 
obeying non-abelian statistics |26], [27] can also be contemplated. Bound states of these kind 
of anyons may be discussed with the help of the semiclassical theory analogous to the one 
considered in this paper. 
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